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Abstract 

This paper consider the possibility of using some quantum tools in decision making 
strategies. In particular, we consider here a dynamical open quantum system helping 
two players, Q\ and Q 2 , to take their decisions in a specific context. We see that, 
within our approach, the final choices of the players do not depend in general on 
their initial mental states, but they are driven essentially by the environment which 
interacts with them. The model proposed here also considers interactions of different 
nature between the two players, and it is simple enough to allow for an analytical 
solution of the equations of motion. 
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I Introduction 


In recent years the scientific literature has seen a growing interest in the possibility of 
using quantum ideas and quantum tools in the description of some aspects of several 
macroscopic systems, systems which, in the common understanding, are usually thought 
to be purely classical. This interest has touched very different fields of science, starting 
with finance, going to ecology, passing through psychology, decision making and so on. 
The literature is now very rich, and it increases almost every day. We just cite here some 
recent books, which cover some of the area mentioned above, but not only. 

In this paper we will propose a dynamical approach to a very simple and well known 
problem in decision making, but in a slightly modified version. Our starting point is what 
was considered in [7], which is a variation on the theme of the prisoners’ dilemma. This is 
just one of the several contributions existing in the literature related to decision making 
processes and to brain dynamics, and it is just one of the contributions suggesting the 
relevance of something quantum in this kind of problems. For instance, Manousakis in 
[8] suggests that the condition describing someone who must still make a choice, could 
be thought as a superposition of suitable states in a particular Hilbert space, whose co¬ 
efficients are related to the probabilities of making a particular choice among the various 
possibilities. He also proposes a time evolution driven by some hyper-simplified hamilto- 
nian. Other effective hamiltonians are used, in similar contexts, by other authors, PHD]. 
Of course these effective hamiltonians, as such, are usually quite ad hoc and can only be 
used to describe some particular aspect of the system under analysis. 

Less dynamically oriented is the paper by Agrawal and Sharda, m, where the au¬ 
thors focus particularly on the probabilistic aspects of the process of decision making. 
Still other contributions are due, for instance, to Vitiello, Khrennikov et ah, [12], and 
to Busemeyer et ah, [13]. In particular, in this last paper the authors confront Markov 
models of human decision-making with other models somehow connected to quantum me¬ 
chanics. A common feature of almost all these papers have to do with the probabilistic 
interpretation of quantum mechanics, where interference effects are quite naturally intro¬ 
duced and described, with respect to what happens using classical ideas, where a similar 
interpretation is not as natural. Other interesting references are given in na. 

Let us now go back to the problem we are interested in here. We adopt almost the 
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same notation as in [7|. We have two players, Gi and G 2 , and each of them can make 
two possible choices, ”0” and ”1”. For instance, 0i means that Qi has chosen ”0”, while 
I 2 means that Gi has chosen ”1”. More compactly, this choice is indicated as O 1 I 2 . 
the same way I 1 I 2 means that both Gi and G2 made the same choice, ”1”. And so ( 
Now, let ns introduce as in [7] four real numbers a, b, c and d satisfying the inequalities 
c> a > d > b. The payoff of Gi is a or c if ^ 2 ’s choice is 0 2 : 0i0 2 corresponds to a payoff 
a, while I 1 O 2 corresponds to c. On the other hand, if Q 2 chooses 1, I 2 , f/i’s payoff is b or 
d: O 1 I 2 corresponds to 6, while lil 2 corresponds to d. The situation is summarized in the 
following table, [7]: 


0 . 


Gl\G2 

0 2 

I 2 

0l 

(a\a) 

(b\c) 

ll 

(c\b) 

(d\d) 


This table represents the point of view of both Gi and For instance, if the two 
players choose 0 (O 1 O 2 ), they both have a payoff a. Analogously, choice I 1 I 2 corresponds 
to the same payoff d for Gi and G2- On the other hand, different choices of the players 
correspond to different payoffs: the choice 0il 2 produces a payoff b for Gi, and c for G2, 
while I 1 O 2 produces a payoff c for Gi, and b for Q 2 ■ The table shows that if Gi chooses 1, 
then he can get the maximum payoff, c (if Q 2 chooses 0), or a small one, d (if Q 2 chooses 
1). Hence this could be the best choice, but it can also have bad consequences. On the 
other hand, if Gi chooses 0, then he can have the minimum payoff, b (if G2 chooses 1), or 
a better one, a (if G2 chooses 0). Hence Gi should make choice 1, hoping that Q 2 makes 
choice 0, so that he gets the maximum payoff c. Of course, with this choice there exists 
also the possibility that he gets d, which is less than a (corresponding to O 1 O 2 ). But, for 
sure, choosing 1, G 1 will not get b, which is the lowest possible payoff. Hence, if what G 1 
really hopes is not to get the lowest payofl§|, he must choose 1. A similar analysis of the 
table from the point of view of G2 suggests that, if also G2 wants to avoid to get the worst 

1 Notice that in any choice X\Y 2 the indices 1 and 2 refer to the players, while X and Y refer to the 
possible choices of the players, 0 or 1. 

2 This is what in the literature is called loss-aversion 
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payoff, he has to choose 1. Then, if Q\ and C/ 2 are rational players , meaning that they are 
both interested not to get b, they should both choose 1 (li 1 2 )- 

In [7j this problem has been considered using quantum techniques: a mental state 
vector belonging to some suitable Hilbert space is associated to each possible choice of 
the players, and it is used to describe the situation. In particular, since we have here only 
four possible choices, the Hilbert space is four dimensional. We will give more details, 
adapted to our aims, in Section II. The dynamics of the vector is deduced by a master 
equation, and the final decision is related to the equilibrium solution of this equation. 

I 11 this paper we consider a similar system from a slightly different point of view, i.e. 
from the point of view of quantum open systems. In our opinion, this choice is more 
realistic, since we consider the possibility that the two players interact with the external 
world, to make up their mind and to take their decisions. Of course, this is different 
from the standard version of the two players game, where there is no interaction at all, 
and this is the reason why we talk of a ”similar system”. In particular, in our case, the 
Hilbert space of the model is richer than the one considered in the existing literature. I 11 
fact, most of the papers considering a quantum approach to decision making deal with 
finite dimensional Hilbert spaces. This is not the case for us: in our settings, while the 
players will be attached to a four-dimensional Hilbert space, the reservoir will not. But, 
rather than being a problem, in our opinion this makes the structure more realistic. In 
fact, the presence of the reservoir mimics well the very many inputs that each player 
normally takes into account while making his choice. This is exactly our interpretation 
of the reservoir: it represents the set of rumors, ideas, suggestions,... coming from the 
external world and reaching the players. For this reason, the dynamics of the players is 
provided by an hamiltonian, written following the rules proposed in [4j, which describes 
not only the two players, but also the reservoir, and, above all, the possible interactions. 
As we have already said, similar ideas have already been used in the literature on decision 
making, see mum for instance. However, in these cases, the hamiltonian is often a very 
simple matrix which, of course, can only be used to describe a particular aspect of the 
model. On the other hand, our hamiltonian contains rather general information on the 
system, and the different ingredients of H can be easily identified. Finally, even if, in our 
knowledge, this is not done in any standard (or quantum) view to the two players game, 
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we will also consider here the possibility of having some interaction between Q\ and Q 2 , 
and we will discuss the consequences of this interaction. In particular, we will consider 
the case in which the two players react in the same way and the case in which they have 
opposite reactions. This will be clarified in the next section. Of course, the presence of 
this interaction between the players, makes our system even more different from the one 
considered in [7], 

The paper is organized as follows: in the following section we propose the model and 
we derive its dynamics. Then we consider the cases in which Q\ and Q 2 do not interact, 
and the cases in which they do, and different possibilities are considered. The analysis of 
the results and our conclusions are discussed in Section III. Finally, to keep the paper self 
contained, we discuss some important facts in quantum mechanics in the Appendix. 

II The model and its dynamics 

In this section we will discuss the details of our model, constructing first the vectors of the 
players, the hamiltonian of the system, and deducing, out of it, the differential equations 
of motion and their solution, with particular interest to its asymptotic (in time) behavior. 

In our game we have two players, Q\ and C/ 2 - Each player could operate two possible 
choices, 0 and 1. Hence we have four different possibilities, which, following [7j, we 
associate here to four different and mutually orthogonal vectors in a four dimensional 
Hilbert space 1-Lg. These vectors are <p 0 ,o, <Pi,o> ‘Pop and <Pi, i- The first vector, <p 0i0 , 
describes the fact that, at t = 0, the two players have both chosen 0 (O 1 O 2 ). Of course, this 
is not a fixed choice, and can change during the time evolution of the system. Analogously, 
ip 0,1 describes the fact that, at t — 0, the first player has chosen 0, while the second has 
chosen 1 (O 1 I 2 ). And so on. = {<Pk,i, k,l = 0,1} is an orthonormal basis for Hg. The 
general mental state vector of the system Sg (i.e. of the two players), for t = 0, is a linear 
combination 

1 

v b = ^ otk,m t h (2-f) 

k,l =0 

where we assume that 1 =0 |o<fc,j| 2 = 1 in order to normalize the total probability. Indeed 
|cco, 0 : is the probability that Sg is, at t = 0, in a state (p 0 ,o, i.e. that both Q 1 and Q 2 have 
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chosen 0. Notice, incidentally, that = $1 {g><h 2 , where <hfc = +x[ k ' > cp{ k \ k = 1,2, 

and where a^i and are related in an obvious way: a 0 ,o — x o >x ^\ 0 : 1,0 = x 'P x i\ 
00,1 = x^xf' 1 and ct 1.1 = x^xf\ We see that the vectors describing Q\ and Q 2 are 
independent, and T is the tensor product of the two. 

The first essential difference with respect to what is done in [7] is now the way in 
which these vectors are constructed: we consider two fermionic operators, see Appendix, 
i.e. two operators b\ and b 2 , satisfying the following canonical anti-commutation rules 
(CAR): 

{bk, b\} = 5k, 1 1, {bk, bi} = 0, (2.2) 

where k, l — 0,1, 1 is the identity operator, and {x,y} = xy + yx. Then we take <p 0 ,o as 
the vacuum of b\ and b 2 : &i<po,o = ^ 2 ^ 0,0 = 0, and construct the other vectors out of it: 

^1,0 — &i¥>0,0, ^0,1 = ^^0,0; ^1,1 — b \ 62^0,0 • 


The explicit expressions of these vectors and operators can be found in many textbooks 
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in quantum mechanics, see pHi] for instance: p>k,i = Vk ® 1 where <p 0 = 

. Then, 


0 


and 


<pi = 



<Pi,o 



®<do 2) 


0 

1 



( 1 ) „ ( 2 ) 
^1,1 = V 5 ! ® 


0 

1 



and so on. The matrix form of the operators bj 


and b\ are also quite simple. For instance, 


bi = 



b 2 = 




and so on. 

Let now hj = b^bj be the number operator of the j-th player: the CAR above imply 
that fiiipk,i — kipk,i and h 2 ifk,i = hpk,h k,l — 0,1. Then, as already stated, the eigenvalues 
of these operators correspond to the choice operated by the two players at t — 0: for 
instance, <pi 5 o corresponds to the choice 1 10 2 , just because one is the eigenvalue of hi and 
zero is the eigenvalue of h 2 . 
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Remark:— One might wonder why, in the description of our model, we use fermionic 
rather than bosonic operators, as we have done in several other applications in recent 
years, [3]. This is easily understood, since the eigenvalues of the fermionic number oper¬ 
ators hj are exactly 0 and 1, which are the only possible choices of the players. On the 
other hand, see [IB] , the eigenvalues of the bosonic number operators are all the natural 
numbers (including 0): too many for us! 

Our main effort now consists in giving a dynamics to the number operators hj , follow¬ 
ing the scheme described in [4], Therefore, what we first need is to introduce a hamiltonian 
H for the system. Then, we will use this hamiltonian to deduce the dynamics of the num¬ 
ber operators as hj(t) := e lHt hje~ lHt , and finally we will compute the mean values of 
these operators on some suitable state which is needed to describe, see below, the status 
of the system at t = 0. The rules needed to write down H are described in [4]. The main 
idea here is that the two players are just part of the full system: in order to take their 
decision, they need to be somehow informed. In fact, it is really the information which 
creates the final decision. Hence, Sg must be open , meaning with this that there must be 
a reservoir 1Z = 1Z\ ® 77-2, interacting with Q\ and C/ 2 , which is responsible for this sort 
of information. The reservoir, compared with Sg, is expected to be a very large system 
since the information is created by several different sources. A possible hamiltonian is 
therefore the following: 

[ h = H 0 + H!, 

i = + (2.3) 

I H, = Y.U V Im (tAW + dk. 

Here ojj and A j are real quantities, and flj(k) are real functions. In analogy with the bj’s, 
we adopt fermionic operators Bj(k) and f?j(fc) to describe the reservoir. They depend on 
j — 1,2 (two different sub-reservoirs for the two players), and on the real variablqj k, and 
they satisfy the rules 

{Bi{k), £,(#} = S itl 5(k - q) 1, {Bi(k),Bj(k)} = 0, (2.4) 

3 In principle we should use a discrete variable to label each element of the reservoirs. However, since 
integrals are quite often easier to be computed than series, as usually done in the literature we consider 
this label to be real. 
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which have to be added to those in (12. 2 j) . Moreover each bj anti-commutes with each 
Bj(k ): {bp B\(k)} = 0 for all j, l and k. Here X ** stands for X or Xb 

II. 1 An interlude: why fermionic operators, and why this hamil- 
tonian? 

It may be useful to recall now that, as discussed in the Appendix, the fermionic operators 
considered above have a very useful characteristic for us: they can be used to construct 
new self-adjoint operators, the number operators h 3 = bjb 3 , which are diagonal in the 
tpk/ s, and whose eigenvalues are exactly zero and one. Which, important to stress, are 
the only two possible choices of our players. As we have already said, this is the core of 
our choice: we have two main possible choices of the players, and these correspond exactly 
to two eigenvalues of very simple matrices. Then, as we have already discussed, a rather 
natural possibility to describe the process of decision making is simply to give a dynamics 
to hj. And our claim is that this dynamics is given (in part) by the hamiltonian (12.31) . 
The full hamiltonian is given below, in (12.51) . 

Let us now concentrate on the meaning of h, beginning with the role of the parameters 
and of the functions. Of course, A j is an interaction parameter, measuring the strength of 
the interaction between Q 3 and 7 Z 3 . If, in particular, Ai = A 2 = 0, then h = H 0 and, since 
[h, hj] = 0, this would imply that the number operators describing the choices of the two 
players stay constant in time. In other words, in this case the original choices of Q\ and Q 2 
are not affected by the time evolution)! Both u 3 and Qj (b) are related to a sort of inertia 
of the system, [4], i.e. to a tendency of a particular part of the system not to change too 
fast its status. For instance, we will see in Section II.1 that Qj(k) is related to the time 
needed by Q 3 to make his choice. In [0] it is also shown, in several concrete applications, 
that the values of uj 3 and klj{k) are related to the magnitude of the oscillations of some 
relevant functions of the model. For this reason, in analogy with classical mechanics, we 
adopt the word inertia in connection with these quantities. 

Let us now explain why we have chosen these particular forms of H 0 , and of Hj. 

4 Please consider that h is not really the full hamiltonian, see (12.511 . In fact, h 3 {t) stays really constant 
in time if [ H , hj] = 0, which is surely true if Ai = A 2 = 0 and if y e x = Hcoop = 0, see (12.511 . 








H 0 is a sum of diagonal operators, describing the free evolution of the operators of 
S = Sg <S> 7Z. In fact, for instance, J2^=i u j^jbj — which is already diagonal 

in terms of the tfk,i- Slightly more complicated, but not particularly different, is the part 
of Ho which refers to the reservoirs: it is also diagonal. Now, suppose that Q\ and Q 2 
are, at t — 0 in a definite state, say (he. O 1 I 1 ), and that the dynamics of the system 
is only given by H 0 . Then, at t > 0, the two players will be still described by <p 0 ,i- no 
change in their decisions. This is coherent with the fact that, as we have discussed before, 
[Ho,hj\ = 0. Stated with different words, we could say that H 0 is the simplest quadratic 
self-adjoint operator in our fermionic operators which commutes with hi and fi 2 - This 
ensures that, in absence of interactions, Q\ and Q 2 do not change idea. 

More interesting is the role of Hi. In order to explain its meaning, we have to recall 
that, see Appendix, bj and B 3 (k) are lowering operators, while their adjoint bj and Zjj (k) 
are raising operators. For instance, if we consider bi<p>i : o, we obtain (p 0 ,o- Then, the action 
of b\ modify the original choice of the players, I 1 O 2 , to the new choice, O 1 O 2 . Similarly, 
since b\ip 0 ,o = the action of brings O 1 O 2 to I 1 O 2 . The operators Bj(k) and Bj(k) 
behave similarly for the reservoir. 

Then, it is clear that Hj describes the interaction between the two components of 
7 Z, Hi and IZ 2 , with the players: bjBj(k) describes the fact that, when the amount of 
information reaching Qj increases (because of B^(k)), Q 3 tends to chose 0 (because of 
bj). On the other hand, B 3 (k)b j describes the fact that Qj tends to chose 1, when the 
amount of information reaching him decreases. Now, recalling that, in our model, what 
the players really want to avoid is getting the smallest payoff 6, and recalling that this is 
achieved by choosing 1, it is natural to interpret the information produced by the reservoir 
as information of bad quality, the more it reaches Qj, the more he moves away from his 
rational choice. 

II.2 Enriching the model 

To make the situation richer and more interesting for us we admit here the possibility 
that the two players also interact among them and we consider two different possible 
interactions, by adding a cooperative and an exchange effects. The full hamiltonian H is 
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therefore 


(2.5) 


H h. hinti 

hint Hex T “1“ Hcoop (^1^2 ^2^1 

where Hex and Hcoop are non negative. In particular, they could be both equal to zero, and 
in this case H = h. In this particular case, Q\ and C/ 2 do not interact with each other. 
On the other hand, if Hex 7^ 0 and Hcoop = 0 then Q\ and Q2 are pushed to make different 
choices because of the terms b\b 2 and 6261, while they act cooperatively if Hex = 0 and 
Hcoop 7^ 0 (because of b\b\ and 6261). Finally, we also allow the possibility of having both 
these contributions, when Hex and Hcoop are simultaneously non zero. 

Before deducing the time evolution of the relevant observables of the system, it is 
interesting to discuss the presence, or the absence, of some integrals of motion for the 
model. In our context, these are (self-adjoint) operators which commute with the ha.mil- 
tonian. In many concrete situations the existence of these kind of operators gives an hint 
on how the hamiltonian should look like, [I], and can be used sometimes to check how 
realistic our model is. In feat, this strategy was previously used to fix the form of H 0 . 
Let us introduce 


2 2 

N =ii N i=ti 

3 =1 3 = 1 




BUk)Bj(k) dk 


( 2 . 6 ) 


with obvious notation. First of all, it is easy to check that [Nj,h] = 0, j = 1,2, so that 


[N, h\ — 0. Moreover, even if 


Nj , p, e 


N, He 


l u 2 


\ u 2 


+ b\b 1 

+ b\bi 


7^ 0, we find that 


= 0, 


so that N commutes also with h + Hex [b\b 2 + b\b\ j. On the other hand, neither Nj nor 

N commute with Hcoop (b\b\ + 62^1 j so that, when Hcoop 7^ 0, IV ceases to be an integral 
of motion. This suggests that the cooperation destroys the integral of motion in (12.bD . 
This is because the cooperative term in H forces Q\ and Q2 to behave in a similar way 
forcing, as a consequence, the mean value of N to change with time. On the other hand, 
if Hcoop = 0, the creation and the annihilation operators in H always compensate their 
actions and for this reason N stays constant in time, even if its different contributions in 
(12.6p have a non trivial time evolution. 
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We can now go back to the analysis of the dynamics of the system. The Heisenberg 
equations of motion X(t) = i[H,X(t)], see Appendix, can be deduced by using the CAR 
(12. 2 p and (12.4p above: 


{ hit) = —iuibi(t) + i\i J R Bi(k,t) dk - ipLexh{t) - i^coo P b\(t), 

hit) = -iu 2 h{t) + i \ 2 J R B 2 (k,t) dk - ifiexhit) + i^coopb\{t), 

Bj(k,t ) = — iQj{k)Bj{k,t) + i\jbj(t), 


j — 1,2. The third equation can be rewritten as 

Bj(k, t ) = + iAj [ bjitJe-^hW-ti) dtl 

Jo 

and, taking Q, j(k) = £ljk, ttj > 0, standard computations produce 


Bj(k, t) dk — [ Bjitye-^dk + in^-bjit). 

Jr ^ ‘‘j 


(2.7) 


( 2 . 8 ) 


We refer to |3] for details of this computation and for a discussion on the physical genesis 
of this approach. If we now replace (12.8ft in the equations (12.7ft for bj(t), we can write 


b{t) = i Ub(t) + i/3{t), (2.9) 

where we have introduced Uj = iojj + / 3j(t ) = J R Bj[k)e~ l ^ kt dk, j = 1,2, and 



( bi{t) A 


/ A Mt) \ 


( iv i 

t^ex 

0 

t^coop 

\ 

b{t) = 

hit) 

b\it) 

, Pit) = 

x 2 p 2 it) 
-Aiftt(t) 

U = 

l^ex 

0 

iv 2 

H'coop 

H'coop 

%V{ 

0 

l^ex 



\ b i(t) / 


--W3 1(0 ) 


\ l^coop 

0 

t^ex 

0*2 

/ 


The solution of (12. 9 p is easily found in a matrix form: 

b{t) = e iUt b{ 0) + i [ e iUit - tl] /3(h) dt u 

Jo 


( 2 . 10 ) 


which is now the starting point for our analysis below. 
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II.3 Q\ and Q 2 do not interact 

This is almost the classical two players game, since they do not interact each other, but 
still both communicate with their environments. As we have discussed before, in this case 
/j, ex — Hcoop — 0. Then U is a diagonal matrix, and e lUt is diagonal as well. Then, from 
(12.1 Op we easily deduce that 

bj(t) = e~ Ujt bj( 0) + i f f3j(t) dti, 

Jo 

j = 1,2. From this equation we can obtain &t(£) and, consequently, the number operator 
fij(t) = bj(t)bj(t). However, what is relevant for us is not really hj(t) itself, but its mean 
value on some suitable state on S. These states are assumed to be tensor products of 
vector states for Sg and states on the reservoir which obey a standard equation, see below. 
More in details, for each operator of the form X s <g) Yji, X s being an operator of Sg and 
Yu an operator of the reservoir, we consider 

(X s ®Y n ) :={*,X S V) uj n (Y n ). 

Here T is the vector introduced in (12.ip . while cun(.) is a state satisfying the following 
standard properties, [4]: 

^ 7 zip'll) = 1, z(Bj(k)) = lu-j z(B^(k)) = 0, ui-ji(B^(k)Bi(q)) = Nj 5jj5{k — q ), (2.11) 

for some constant Nj. Also, 0 Jn{Bj(k)Bi(q)) = 0, for all j and l. These formulas for 
reflect for the reservoir expressions similar to those for Sg. Then 

Ujit) = ( hj(t )) = e- 27rA ? /Q A||^||2 + Nj _ e -2 nfyfljt') _ (2.12) 

What is interesting here is that, if X 3 7 ^ 0, 

rij( 00 ) := lim rij(t) = Nj (2.13) 

t—t OO 

does not depend on the original state of mind of the two players, but only on what the 
reservoir suggests. In fact, independently of the vector T describing probabilistically, at 
t = 0 , the choices of both Q\ and if 1Z.\ (the part of the reservoir interacting with Qi) 


12 






has in (12.111) Ni = 0, then after a sufficiently long time, 0 will be exactly QQs choice. On 
the other hand, if Ny = 1, then Qy will eventually choose 1. A similar conclusion can be 
deduced for Q 2 . Therefore, when Qy and Q 2 do not interact, their choices are only dictated 
by their environments. This conclusion looks quite reasonable, in the present context. 

Remarks:— (1) Notice that, if A j = 0, formula (I2.12p reduces to rij(t ) = ||fcjT|| 2 = 
rij( 0), V t. This is not surprising, since reflects what was already deduced before in absence 
of interactions of any kind. In this case, in fact, we have seen that the initial state of 
mind is what really matters for the final decision, since there is no time evolution of the 
operator hj at all. 

(2) More in general, formula H2.12j) suggests the introduction of a sort of characteristic 
time for Qj, Tj = The more t approaches Tj, the bigger the influence of IZj on Qj is. 
In particular, if \ 3 —> 0, Tj diverges. Hence we recover our previous conclusions: Qj is not 
influenced at all by 7 Zj, even after a long time. A similar behavior is deduced also when 
Qj increases: the larger its value, the larger the value of Tj. In other words, for large Qj 
the influence of the environment is effective only after a sufficiently long interval. This 
is not very different from what we have deduced in other systems, [3], where analogous 
parameters of the hamiltonian measure the inertia of that particular part of the system. 
Of course, r 3 can be considered as a sort of decision time. 

(3) Since the rational choice of both players is 1, (12.131) shows that rationality really 
belongs to IZj , rather than to Qj: in our version of the game, Qj does not need to be 
rational, at least if their reservoirs behave rationally! 

II.4 The effect of exchange interaction 

In the following we will fix // coop = 0, allowing /j, ex to be different from zero. In particular, 
from now on, for concreteness’ sake we will work fixing the following values of the other 
parameters in the hamiltonian: ojy = 1, co 2 = 2, Xy = A 2 = 0.5, Qx = Q 2 = 0.1. This 
choice is meant to have almost identical players and reservoirs. As it is clear, the only 
difference between Qy and Q 2 is played here by the values of uq and ui 2 . 

After few computations, calling V(t) = e lUt and 14 ,i{t) its (k, l) -matrix element, we 
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deduce that 


ni(t) = |y 1 ,i(t)| 2 ||6i^|| 2 + |y 1 , 2 (i)| 2 ||^|| 2 + 


27r / dti 

Jo 


A? 


A| 


ni IU,i(t - ii)| 2 iVi + IU,2(i - ii)| 2 JV 2 


(2.14) 


and 


n 2 (t) = |^ 2 ,i(t)| 2 |MI| 2 + |V 2 , 2 W| 2 ||^|| 2 + 


27 t f dt\ 
Jo 


- h)] 2 ^ + ^|y 2)2 (t - ti)| 2 iv 2 


(2.15) 


To begin with, we consider three different choices for /j, ex : (a). fi ex = 0.01, (b). fi ex = 0.05 
and (c). n ex — 0.1. In all these cases it is possible to check that both VAi (t) and V\ i2 (t) 
converge to zero when t diverges. On the other hand, neither f* dt\\ VX,i(t — ti)\ 2 nor 
dti\Vi t2 (t — ti)\ 2 converge to zero. All these computations can be performed analytically 
and the explicit result, in case (a), is the following: 

ni(oo) ~ 0.99997Ah + 0.000017V 2 , 
n 2 (oo) ~ O.OOOOllVi + 0.99997-/V 2 . 


(2.16) 


How we can see, these are symmetrical, and not very different from the result in (|2.13[) : 
the two players modify their decision with respect to when fi ex = 0, but just a little bit ! 
This is because fi ex is too small. In fact, let us consider the case (b) above, /j, ex = 0.05. 
In this case, repeating the same computations, we conclude that 


ni(oo) ~ 0.99997iVi + 0.00251iV 2 , 
n 2 (oo) ~ 0.00251#! + 0.999977V 2 , 


(2.17) 


which shows that the mixing between N\ and N 2 increases a little bit. And, in fact, this 
mixing increases even more in case (c), when fi ex = 0.1: we get 


n 1 ( 00 ) ~ 0.99039iVi + 0.00958iV 2 , 
n 2 (oo) ~ 0.00958^ + 0.99039# 2 . 


(2.18) 


To clarify further the role of the exchange hamiltonian, we now consider much higher 
values of fi ex , keeping again /i coop = 0. Hence we take: (d). / i ex = 10 and (e). fx ex = 100. 
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In the first case, fi ex = 10, we find 


n 1 (oo) ~ 0.596277Vi + 0.403707V 2 , 
n 2 (oo) ~ 0.40370AA + 0.596277V 2 , 


(2.19) 


while in case (e), fr ex = 100, we obtain 


ra^oo) ~ 0 . 50154 ^ + 0.498467V 2 , 
n 2 (oo) ~ 0.498467V! + 0.501547V 2 . 


( 2 . 20 ) 


We believe that, for p, ex \i coop = 0, the two players reach eventually a common choice 
which should be n^oo) = ro 2 (oo) = 1(7^ + 7V 2 ): perfect mixing! Once again, then, the 
decisions of Q\ and Q 2 are driven by the reservoirs but, in this case, the stronger the 
interaction between Qi and (? 2 , the more TZ\ and 1Z 2 affect in a symmetric way the two 
players. 

II.5 The effect of cooperative interaction 

We now consider the case in which only the cooperative part in the hamiltonian is switched 
on, /i coop 7^ 0, while the exchange contribution is turned off, fi ex = 0. As before, we 
will consider, for the same reasons, the following values of the other parameters in the 
hamiltonian: — 1, 0 J 2 — 2, Ai = A 2 = 0.5, Ox = 0 2 = 0.1, and then we will put (a). 

t^coop 0.01, (b). t^coop 0.05 and (c). ia COO p 0.1. 

In this case we deduce that 



J 0 


ni(t) = |I/i, 1 (t)| 2 ||5 1 4/|| 2 + |Vi i4 (f)| 2 (1 - ||& 2 *|| 2 )+ 
[dh 4l|V,,i (t - + A| v L2it - tl )| 2 (l - JV 2 ) 


( 2 . 21 ) 


and 


n 2 (t) = |V2 i2 (t)| 2 ||5 2 'h|| 2 + \V 2 , 3 (t)\ 2 (l - ||6 1 'h|| 2 ) + 



( 2 . 22 ) 
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Again it is possible to check that all the functions 14 ,i(t) above converge to zero when 
t diverges. On the other hand, f {) dt\\Vkj(t — ti )| 2 admits a non zero limiting value for 
t —> 00 . The results are the following: in case (a), /i COO p = 0.01 and p ex = 0, we have 


J ni(oo) ~ 0.99997Ah +0.00001(1 -Ah), 
| n 2 (oo) ~ 0.00001(1 - Ni) + 0.99997N 2 . 

In case (b), p. coop = 0.05 and p, ex = 0, we have 

f m ( 00 ) ~ 0.9996677! + 0.00028(1 -Ah), 
{ n 2 ( 00 ) ~ 0.00028(1 - Ah) + 0.99966Ah, 

while in case (c), ^ coop = 0.1 and fi ex = 0, we get 

f ni(oo) ~ 0.99887A^i + 0.00115(1 — N 2 ), 
\ n 2 ( 00 ) ~ 0.00115(1 - Ni) + 0.99887^2. 


(2.23) 


(2.24) 


(2.25) 


Again we observe that the higher the value of ^ COO p, the higher the mixing between the 
effects of the two sub-reservoirs. Hence we are led to formulate a similar conclusion as we 
did in the previous situation, and we expect that, for + coop /j, ex = 0, the two players 
arrive to an asymptotic (in time) choice which is the following: ni(oo) = \(Ni + (1 — N 2 )), 
n 2 ( 00 ) = |((1 - Ah) + Ah). 


II.6 Full hamiltonian 

In this last part, we consider together the effects of the exchange and of the cooperative 

hamiltonians, still keeping unchanged the values u\ — 1, oj 2 — 2, Ai = A 2 = 0.5, Oi = 

fl 2 = 0.1. Now both fa ex and ^ coop will be taken different from zero. In particular, we will 
consider the situation in which fi ex and /i COO p are significantly different from each other 
(which we don’t expect is particularly different from what we did before), and the case 
in which they are similar. More in details, these will be our choices of parameters: Case 
(a). p<ex 0.01 and ia COO p 100, (b). f-i ex 0.01 and ia COO p 1, (c). f-i ex t-^coop 0.5, 

(h). fa ex 1 and P’coop 0.01 and (e). ia& x 100 and Pcoop 0.01. 

In this case we deduce that 

^(t) = |F lil (t)| 2 || 6 1 4>|| 2 + |H li 2 (t)| 2 || 6 2 4>|| 2 + |H 1 ) 3 (t)| 2 (l-|| 6 1 4 /|| 2 ) + |H M (t)| 2 (l-l| 6 2 4 >|| 2 )+ 
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(2.26) 


and 


+27r 



J A 2 


[\v 1A (t - ti)| 2 iVi + \V 1>3 (t - ti)| 2 (l - iVi)] + 


+ 27T 



[I v 1)2 (t - h) \ 2 N 2 + | V 1A (t - t\) | 2 (1 - N 2 )] 


^(t) = |^ 1 (t)| 2 ||6 1 ^|| 2 + |y 2i2 (t)| 2 ||6 2 ^|| 2 + |y 2 , 3 (t)| 2 (l-||^|| 2 ) + |y 2i 4(t)| 2 (l-||6 2 ^|| 2 )+ 


+27r f [| V2,i(t ti)| 2 iVi + | V 2j3 (t — ti)| 2 (l — iVi)] + 

Vo “i 

+ 2tt [ dt^ [|F 2 , 2 (t-t 1 )| 2 7V 2 + |F 2)4 (t-t 1 )| 2 (l-iV 2 )] . 
Vo “2 


(2.27) 


The following are the results we have deduced in the five cases listed above. We have: 
Case (a), /i ex = 0.01 and /i coop = 100: 


ni(oo) ~ 0.50317W + 0.49682(1 - N 2 ), 
n 2 (oo) ~ 0.49682(1 - N x ) + 0.503171V 2 . 


(2.28) 


Case (b), fi ex = 0.01 and /x coop = 1: 


n 4 (oo) ~ 0.91914W + 0.08075(1 - N 2 ), 
n 2 (oo) ~ 0.08075(1 - AT,) + 0.91914iV 2 . 


(2.29) 


Case (c), p e3; = 0.5 and /i coop = 0.5: 


' m(oo) ~ 0.85428W + 0.00626(1 - ,¥,) + 0 . 119747 V 2 + 0.01917(1 - N 2 ) = 
= 0.84802W + 0.10057iV 2 + 0.02543, 
n 2 (oo) ~ 0.11974W + 0.01917(1 - 1V 4 ) + 0.854281V 2 + 0.00626(1 - N 2 ) = 
= 0.10057W + 0.84802 W + 0.02543. 


Notice that, in these equations, the first form has been explicitly written simply because, in 
this way, the different contributions arising from (12.2611 and (12.271) can be easily identified. 
Case (d), /i ex = 1 and /i coop = 0.01: 


ni(oo) ~ 0.99210W + 0.007957V 2 , 
n 2 (oo) ~ 0.00795W + 0.99210iV 2 , 


(2.31) 


17 








and, finally, Case (e), fx ex = 100 and fi COO p = 0.01: 

I m(oo)~ 0.50308^ + 0.49692^, 

[ n 2 (oo) ~ 0.49692W + 0.50308iV 2 . 

We postpone our detailed analysis of these and of the previous results to the next section. 
Here we just want to add that, contrarily to what we have seen in Section 11.1, in this 
more general case we expect that the characteristic time depends also on fj, ex and /i coop , 
so that these parameters are expected to contribute to the decision time. 


Ill Analysis of the results and conclusions 

The first clear output of our analysis suggests that, when Q 1 and Q 2 do not directly 
interact, it is really the environment which produces their decisions. Hence the rationality 
of the players is strongly linked to the nature of the reservoirs: if both reservoirs have 
Nj — 1, j — 1,2, then rij( oo) = 1, and the two players make the most rational choice 
according to the loss aversion rule. More interesting is the situation when we allow some 
interaction between Q\ and Q 2 . In particular our results show that, when at least one 
between fi ex or fi coop is different from zero, and small, the value of rij(oo ) is essentially 
decided again by the j-th part of the reservoir. However, when the numerical values of 
one of the two parameters increase, then some mixing is possible. For instance, we see 
that when fi ex = 100 and /a coop = 0, ni(oo) ~ 0.50154W + 0.498461V 2 and n 2 (oo) ~ 
0.49846W + 0.501541V 2 . This means that, even if the two components of the reservoir do 
not mutually interact, the existence of a direct interaction between Q\ and Q 2 mixes the 
cards : the final decision of each player is not only related to the value of his own part 
of reservoir (i.e. to N\ or to N 2 ), but it is a mixture of the two, and, at least for this 
high value of fi ex , in this mixture Ni and N 2 have almost the same weights for Qi and Q 2 . 
A similar behavior is observed also when g, ex = 0 while n coop increases: again we have a 
stronger and stronger mixing of the effects of TZi and 1Z 2 for /i COO p increasing. However as 
we see from (12T23H - (12725H . Nx mixes with 1 - N 2 (rather than with N 2 ) and N 2 with 1-Ni 
(rather than with Ni). Hence, this contribution in the hamiltonian, behaves differently 
from the other one, and this is natural, due to the different kind of the interactions. When 
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we consider both contributions in h int , the two effects come together and we see this in 
formulas (I2.28p - (I2.32I) . From these formulas we also see that in the extreme situations 
(when n ex is much smaller or much larger than fi COO p), not unexpectedly the two final 
decisions of Q\ and Q 2 are similar to the previous cases (i.e. to the cases in which one of 
the /Ts was zero). On the other hand, when / i ex = fi coop , the two effects are both clearly 
visible, see formula (I2.80p . 

In order to compare these results with those in the Introduction, we begin with a very 
evident fact: the initial state of mind of the players plays absolutely no role in the final 
decision, except when there is no interaction at all. No matter which was their status 
at t — 0, its effect simply disappears when t increases. This is clearly a measure of the 
fact that our model is not really the two-player game proposed in [TJ , as we have already 
stressed before, but a slightly different version of that. 

Let us now consider four different cases, depending on the values of Nj of 7 Zj. Case 
(I): Ni = N 2 = 0; Case (II): N\ — 0 and N 2 = 1; Case (III): Ni — 1 and N 2 = 0; Case 
(IV): Ni = N 2 = 1. From the formulas of Section II we deduce the following: 

1. the only way in which both Q\ and Q 2 choose 1 if when N\ — N 2 = 1, but not for 
all values of / i ex and /2 COO p. For instance, apparently this is not so when [i coop <C fi ex . 
However, when this happens, 72,1(00) and 722(00) still coincide. 

2. on exactly the opposite side, when N\ = N 2 = 0 the values of 721(00) and 722(00) stay 
always very low, except again when fi coop <C n ex . Even now, when this happens, 
721(00) and 722(00) still coincide. The larger /i COO p with respect to fa ex , the bigger 
the value 721(00) = n 2 (oo) which approaches asymptotically, as our computation 
suggests, the value |. 

3. when N\ — 0 and N 2 — 1 in most of the cases considered here 721(00) stays close 
to 0 while 722(00) is close to 1. However, when \i ex <C /i CO op, again our numerical 
results suggest that 721(00) = 72 2 (oo) ~ Specular (and similar) conclusions can 
be deduced when N x — 1 and N 2 = 0. 

4. while there is apparently no other way to get 721(00) = 722(00) = 1 than having 
A^i = IV2 = 1, there exist several possibilities to have 721(00) = 72 2 (oo). Therefore, 
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in a slightly modified version of the game in which we look for equal decisions (not 
necessarily equal to 1), we have plenty of possibilities in which this happens. 

Remarks:— (1) A different possibility, which may be closer to the usual interpretation 
of what is quantum in decision making, is to look at the nj( oo) we have deduced before in 
a probabilistic way. For instance, rather than looking at rq(oo) as the real decision taken 
by Qj, we could consider it as a sort of probability that Q 3 chooses 0 or 1. Then, instead 
of looking to square modula of the coefficients of the vectors in T, we directly look at 
rij(oo). But this does not fit well with our general interpretation, see [4], and we will not 
insist on it here. 

(2) It should probably be stressed that the payoffs a, 6, c and d do not enter explicitly 
in the definition of the hamiltonian, at least in the model considered here. In fact, 
we are interested here in the possibility that Qj make the rational choice for a fixed 
choice of parameters satisfying c > a > d > b, whatever this choice is. Changing their 
values, but maintaining these inequalities, we don’t affect the players’ behavior, of course. 
Nevertheless, it could be interesting to look for some different model in which the role 
of the payoffs is evident in the hamiltonian of the system itself or directly in the state 
describing the system at t — 0. The (probably) easiest way to include the payoffs directly 
in the hamiltonian is to assume, for instance, that fi ex and /r coop depend explicitly on a, 
b , c and d. For instance, if this dependence is such that fi ex > n C oo P , then the effect of the 
exchange interaction would be stronger than that of the cooperative term in hi nt . More 
sophisticated dependencies could be considered, like for instance some nonlinear extra 
term in H , depending on the payoffs. But this would make very hard, if not impossible, 
to get an exact analytical solution, and perturbative expansions should be possibly used. 

This is probably just the beginning of the story: there are still several possible aspects 
to be considered. First of all, we have considered here just a particular choice of the 
many parameters of H. A natural question is what changes when these parameters, and 
in particular those which reflect the nature of the players, are fixed in a different way. 
For instance, in view of the meaning of the uij s we have deduced for other systems, [4], 
we could expect a larger inertia of, say, Q\ with respect to Q 2 if 0 J 2 '- G 1 changes his 
original idea slowly, when compared to However, from the point of view of ni(oo) and 
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77,2(00), we don’t expect this will change much our conclusions, but at most some (minor) 
details, like the decision time. Moreover, the hamiltonian we have considered here is just 
one among all the possible choices. Indeed, in our opinion, it is a rather natural choice and, 
when compared with other possibilities, allows a more natural interpretation. Still, one 
could look for other possibilities, and for instance one could try to add non linearities in 
the model. However, in this case, numerical techniques should most probably be adopted. 
Another interesting aspect is the following: is there any other problem in decision making 
theory in which the method proposed here could be applied? We believe this is very 
plausible. These are some of the aspects we plan to consider in a close future. 
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Appendix: Few results on the number representation 

To keep the paper self-contained, we discuss here few important facts in quantum mechan¬ 
ics and in the so-called number representation. More details can be found, for instance, 
in [HU 15], as well as in [4]. 

Let % be an Hilbert space, and B{H) the set of all the (bounded) operators on H. 
Let S be our physical system, and 21 the set of all the operators useful for a complete 
description of S, which includes the observables of S. For simplicity, it is convenient (but 
not really necessary) to assume that 21 coincides with B{H) itself. The description of 
the time evolution of S is related to a self-adjoint operator H = H t which is called the 
Hamiltonian of S, and which in standard quantum mechanics represents the energy of S. 
In this paper we have adopted the so-called Heisenberg representation, in which the time 
evolution of an observable X 6 21 is given by 

X{t) = exp(zL/f)Xexp(— iHt), (A.l) 
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or, equivalently, by the solution of the differential equation 

J - = iexp(iHt)[H, X] exp(— iHt) = i[H,X(t)\, (A.2) 

where [A, B] := AB — BA is the commutator between A and B. The time evolution 
defined in this way is a one-parameter group of automorphisms of 21. 

An operator Z e 21 is a constant of motion if it commutes with H. Indeed, in this 
case, equation (1A.2I) implies that Z(t) = 0, so that Z(t) — Z for all t. 

In some previous applications, |4], a special role was played by the so-called canon¬ 
ical commutation relations. Here, these are replaced by the so-called canonical anti¬ 
commutation relations (CAR): we say that a set of operators {a^, a\,£ = 1,2 ,..., Lj 
satisfy the CAR if the conditions 

{a e , a f n } = 6 £n 1, { a £ , a n } = {a\, a*} = 0 (A.3) 

hold true for all £, n — 1,2,... ,L. Here, 1 is the identity operator and {x, y} := xy + yx 
is the anticommutator of x and y. These operators, which are widely analyzed in any 
textbook about quantum mechanics (see, for instance, mm are those which are used 
to describe L different modes of fermions. From these operators we can construct hi = a\ai 
and N = Y2i= i which are both self-adjoint. In particular, hi is the number operator for 
the £-t\i mode, while N is the number operator of S. Compared with bosonic operators, 
the operators introduced here satisfy a very important feature: if we try to square them 
(or to rise to higher powers), we simply get zero: for instance, from (1A.3|) . we have a“j = 0. 
This is related to the fact that fermions satisfy the Fermi exclusion principle |15j . 

The Hilbert space of our system is constructed as follows: we introduce the vacuum 
of the theory, that is a vector (po which is annihilated by all the operators ap. aitpo = 0 
for all £ = 1, 2,..., L. Such a non zero vector surely exists. Then we act on ip 0 with the 
operators a\ (but not with higher powers, since these powers are simply zero!): 

< -Pni,n 2 ,... 1 nL '■= (4)^ (<4)" 2 ' ' ' (4)^0, (A.4) 

Hi = 0,1 for all £. These vectors form an orthonormal set and are eigenstates of both hi 
and N: hpp ni , n ^ nL = npp ni , n ^ nL and N(p nun2 ,..., nL = Nip n ^ n2 ^ nL , where N = Y^t=\ n t- 
Moreover, using the CAR, we deduce that 

Tli (®£^ni,n2,...,n^) (jl£ 1) (®f^ni,n2,...,n^) 
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and 


\^ l £^Pni,n2,...,ni J J (j^£ ~f" 1) {fll VYi ,n2,...,rii ) j 

for all l. Then an and a\ are called the annihilation and the creation operators. Notice 
that, in some sense, a\ is also an annihilation operator since, acting on a state with 
ri£ = 1, we destroy that state. 

The Hilbert space FL is obtained by taking the linear span of all these vectors. Of 
course, FL has a finite dimension. In particular, for just one mode of fermions, dim(FL ) = 2. 
This also implies that, contrarily to what happens for bosons, all the fermionic operators 
are bounded. 

The vector (p ni ,n 2 ,...,n L in (1A.4H defines a vector (or number) state over the algebra 21 
as 

00 = (v ni,n2,...,n L i Xtp n \,n2,...,nL ), (A.5) 

where (, ) is the scalar product in FL. As we have discussed in |4], these states are useful 
to project from quantum to classical dynamics and to fix the initial conditions of the 
considered system. 
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